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Abstract 

Many limits are known for hypergeometric orthogonal polynomials 
that occur in the Askey scheme. We show how asymptotic represen- 
tations can be derived by using the generating functions of the poly- 
nomials. For example, we discuss the asymptotic representation of the 
Meixner-Pollaczek, Jacobi, Meixner, and Krawtchouk polynomials in 
terms of Laguerre polynomials. 

1 Introduction 

It is well known that the Hermite polynomials play a crucial role in certain 
limits of the classical orthogonal polynomials. For example, the ultraspher- 
ical (Gegenbauer) polynomials C2(x), which are defined by the generating 
function 

oo 

(1 - 2xw + w 2 )-"< = C2(x)w n , -1<x<1, \w\<l, (1) 

n=0 

have the well-known limit 

lim r nl2 Cl{x/^j) = -H n (x). (2) 
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For the Laguerre polynomials, which are defined by the generating func- 
tion 

oo 

{l- w )- a - l e- wx ^ l - w ^ = Y,Ln(x)w n , \w\<l, (3) 

n=0 

a,x G <D, a similar results reads 




These limits give insight in the location of the zeros for large values of the 
limit parameter, and the asymptotic relation with the Hermite polynomials 
if the parameters 7 and a become large and x is properly scaled. 

Many methods are available to prove these and other limits. In this paper 
we concentrate on asymptotic relations between the polynomials, from which 
the limits may follow as special cases. 

In Q many relations are given for hyper geometric orthogonal polyno- 
mials and their q— analogues, including limit relations between many poly- 
nomials. In Figure 1 we show examples for which limit relations between 
neighboring polynomials are available, but many other limit relations are 
mentioned in ||, §, |] and @. 

In m and we have given several asymptotic relations between 
polynomials and Hermite polynomials. In these first papers we consid- 
ered Gegenbauer polynomials, Laguerre polynomials, Jacobi polynomials, 
Tricomi-Carlitz polynomials, generalized Bernoulli polynomials, generalized 
Euler polynomials, generalized Bessel polynomials and Buchholz polynomi- 
als. 

The method for all these cases is the same and we observe that the 
method also works for polynomials outside the class of hypergeometric poly- 
nomials, such as Bernoulli and Euler polynomials. 

Our method is different from the one described in Q and |Q, where also 
more terms in the limit relation are constructed in order to obtain more 
insight in the limiting process. In these papers expansions of the form 

00 

P n (x;X) = J2 R k(x;n)\- k 
k=0 

are considered, which generalizes the limit relation 

lim P n (x; A) = R (x;n), 

A^oo 

and which gives deeper information on the limiting process. In [Q] a method 
for the recursive computation of the coefficients Rk(x;n) is designed. 



2 



In jjjj similar methods are used, now in particular for limits between clas- 
sical discrete (Charlier, Meixner, Krawtchouk, Hahn) to classical continuous 
(Jacobi, Laguerre, Hermite) orthogonal polynomials. 

Askey Scheme of Hypergeometric 
Orthogonal Polynomials 



4^3 



3 F 2 



2 F l 



Wilson 
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Continuous 
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Continuous 
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Laguerre 
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Racah 
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Dual Hahn 

n, x, y, S, N 




Krawtchouk 

n, x, p, N 



2 F 



Charlier 



n, x, a 



Hermite 



II, X 




Figure 1. The Askey scheme for hypergeometric orthogonal polynomials, 
with indicated limit relations between the polynomials. 

In current research we investigate if other limits in the Askey scheme 
can be replaced by asymptotic results. Until now we verified all limits from 
the third level to the fourth (Laguerre and Charlier) and the fifth level 
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(Hermite). Several limits are new, and all results have full asymptotic ex- 
pansions. 



2 Asymptotic representations 

Starting point in our method is a generating series 



F(x,w) = J2Pn(x)w n , (5) 

n=0 

F is a given function, which is analytic with respect to w at w = 0, and p n 
is independent of w. 

The relation (|5|) gives for p n the Cauchy-type integral 

Pn( x ) = ^— I F(x,w" 



2ni Jc w n+1 ' 

where C is a circle around the origin inside the domain where F is analytic 
(as a function of w). 
We write 

F{x,w) =e Aw ~ Bw2 f(x,w), 
where A and B do not depend on w. This gives 

Pn(x) = ^-I ^-^f( x , w )-%. (6) 

Because / is also analytic (as a function of w), we can expand 

oo 

f(x,w) = e- Aw+Bw2 F(x,w) = £) <%«;*, (7) 

k=0 

that is, 

1 



/(x, to) = 1 + \pi(x) - A]w + 



p 2 {x) -Api{x) +B + ^A 2 



w" + 



if we assume that po{x) = 1 (which implies Co = 1). 

We substitute (0) in (|6|). The Hermite polynomials have the generating 
function 



n=0 
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which gives the Cauchy-type integral 

H n (x) = ^- ( e 2xz ~ z2 z~ n ~ l dz, (8) 
2-ki Jc 

where C is a circle around the origin and the integration is in positive direc- 
tion. The result is the finite expansion 

Pn{x)=Z — , Z = VB, e=^=, 0) 

because terms with > n do not contribute in the integral in (|6|). 

In order to obtain an asymptotic property of (^) we take A and B such 
that c\ = C2 = 0. This happens if we take 

A = p 1 (x), B = -p\{x) -p 2 (x). 

As we will show, the asymptotic property follows from the behavior of 
the coefficients Ck if we take a parameter of the polynomial pk(x) large. We 
use the following lemma, and explain what happens by considering a few 
examples. 

Lemma 2.1 Let 4>{w) be analytic at w = 0, with Maclaurin expansion of 
the form 

4>{w) = fJ.w n (a + aiw + a 2 w 2 + ...), 

where n is a positive integer and are complex numbers that do not depend 
on the complex number p, a$ ^ 0. Let Ck denote the coefficients of the power 
series of f(w) = e^ w \ that is, 

oo 
k=0 

Then cq = 1, cj^ = 0, k = 1, 2 . . . , n — 1 and 

c k = o(\^ k ^), ^oo. 
Proof. The proof follows from expanding 

k=0 k=0 
00 k kn 

— - — (a + aiw + a2W + . . .) , 
k=o fc - 

and comparing equal powers of w. ■ 
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2.1 Ultraspherical polynomials 

The generating function is given in ([[]), and we obtain 

A = Clix) = 2x 7 , B= l - [C~{(x)\ 2 - C2(x) = 7 (1 - 2x 2 ). 
The expansion reads 

^-'l?^' (10) 

where z = VtP ~~ 2x 2 ), £ = xy/z. We have 

2 

c = 1, ci = c 2 = 0, c 3 = -7x(4x 2 - 3). 

Higher coefficients follow from a recursion relation. 
The function f(x,w) of @ has the form 

/(x, w) = e ^ x,w \ 4>(x, w) = 7^ 3 (a + a\w + a 2 w 2 + ...)• 

By using Lemma 2.1 and £ = 0(y/^y) we conclude that the sequence {(f>k} 
with = Cfe/ ' z k H n ^k(0 has the following asymptotic property: 



( 7 n/2+|*/3j-*) , /c = 0, 1, 2, 



This explains the asymptotic nature of the representation in (10) for large 
values of 7, with x and n fixed. 

To verify the limit given in (Q), we first write x in terms of £: x = 
£/ + 2^ 2 - With this value of x we can verify that Ck/z k = o(l), 7 — > 00, 
and in fact we have the limit 

lim - — — tkCI I — I = —H n (x). 

7-00 ( 7 + 2x 2 )"/2 » \y 7 + 2x 2 / n! nV ; 

2.2 Laguerre polynomials 

We take as generating function (see (g)) 

00 

F(s,u>) = (l + u,)-«- 1 e W(i+«') = J2(-l) n L%(x)w n . 

n=0 
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We have A = x— a — 1, B = x — \{a + 1), and we obtain 



where z = y^x — (a + l)/2, £ = (x — a — l)/(2z). The first coefficients are 

1 . 

c = 1, ci = c 2 = 0, c 3 = ~{3x - a- 1). 

Higher coefficients follow from a recursion relation. The representation in 
has an asymptotic character for large values of \a\ + \x\ . It is not difficult 
to verify that the limit given in (Q) follows from (pd|). 

3 Expansions in terms Laguerre polynomials 

We give examples on how to use Laguerre polynomials for approximating 
other polynomials. The method for the Hermite polynomials demonstrated 
in the previous section can be used in a similar way. 

Lemma 3.1 Let the polynomials p n {x) be defined by the generating func- 
tion 

oo 

F(x,w) = ^2p n (x)w n , 

n=0 

where F(x,w) is analytic in w = and F(x,0) = 1. Let 

f{x,w) = e -Aw/{B W -i)^ _ Bw) c+1 F{x,w), 
and let the coefficients Ck(x) be defined by the expansion 

oo 

f{x, w) = Y^ c k (x)w k , c = 1, (12) 

fc=0 

where A,B and C do not depend on w. Then p n {x) can be represented as 
the finite sum 

Pn(x)=B^± C -0±L^ k{O , t = ^ (13) 
where L®(x) are the Laguerre polynomials. 
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Proof. The polynomials p n (x) can be written as 

Pn{x) = _L / e Aw /^Hl-Bw)- c - 1 f(x,w)- d ' V 



2-ki Jc ' w n+1 ' 

where C is a circle around the origin in the domain where F(x,w) is analytic 
(as a function of w). By substituting the expansion of f(x, w) and using the 
generating function (||) of the Laguerre polynomials the proof follows. ■ 

This time, A, B and C can be chosen such that c\ = 0, C2 = 0, C3 = 0. 
These coefficients are given by 

Cl = Pl - BC - B + A, 

c 2 = P2-piBC- Pl B+ Pl A-ABC + ^(B 2 C 2 + B 2 C + A 2 ), 

c 3 = P3-p2BC-p2B+p 2 A- Pl ABC-l(B 3 C 3 + B 3 C + A 3 ) + 

o 

^( Pl B 2 C 2 + v\B 2 C + V \A 2 + AB 2 C 2 - AB 2 C - BA 2 C + A 2 B). 

We see that the equations c\ = 0, c 2 = 0, C3 = for solving for A, B and C 
are nonlinear. However, solving c\ = 0, C2 = for A and C gives 

A = B(C+ 1) -pi, C= ^ ' 

and with these values C3 becomes 

cs=P3- P2P\ + ^(Pi^ 2 +Pi + 2 Pi B ~ 4p 2 B), 

and C3 = is a quadratic equation for B. 

As follows from the above representation of C, this quantity will depend 
on x. This gives an expansion for p n (x) in terms of Laguerre polynomials 

(0 with the order depending on x. When studying properties of p n (x) 
(for example investigating the zeros) this may not be very desirable. In that 
case we can always take C = a (not depending on x), and concentrate on 
two equations c\ = 0, C2 = for solving A and B. This gives 

A=^p 2 -(a + l)(2p2-p 2 1 ), B = P -^±^. (14) 

The order a may be chosen conveniently, without requiring C3 = 0. 

For large values of certain parameters in p n {x) expansion (13) may have 
an asymptotic property when taking c\ = C2 = C3 = 0, but also when only 
ci = or ci = C2 = 0. In the following section we give four examples, for 
one level of the Askey scheme, namely for the Meixner-Pollaczek, Jacobi, 
Meixner, and Krawtchouk polynomials. 
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4 Expanding Meixner-Pollaczek into Laguerre po- 
lynomials 



For the Meixner-Pollaczek polynomials we have the generating function: 

-X+ix / - M \ —X—ix 

n=0 

The expansion for the Meixner-Pollaczek polynomials reads 



—X-\-ix —X—ix 

F(x,w) = (l-e^w) "(l-e-^w) " = ^P| A) (x;^K. (15) 



Pi x \x;<f ) ) = Y / B n ~ k c k L^ k (0, S = A/B, (16) 

k=0 



where the coefficients c k follow from (|12| ) with F(x,w) given in (|15|). 

We write x + iX = re ie ,9 G [0, tt],t > 0, and consider r — ► oo; the 
asymptotic results hold uniformly with respect to 9. 

4.1 One free parameter 

First we consider a simple case by taking B = 1 and C = a, and solve ci = 
for A. This gives 

.A = a + 1 — 2A cos <p — 2x sin 0. 
The first coefficients are given by 

Co = 1, ci = 0, C2 = x sin 2</> + A cos 2</> — 2(x sin + A cos 4>) + -a, 
and the remaining ones can be obtained from the recursion 

(k + l)c k+ i = 2(1 + cos (f>)kck + 

[a + 1 — 2A + 4(cos (j) — 1) (A cos ^ + x sin 0) + 
2(1- A:)(l + 2cos0)]c fc _i + (17) 
[4A + 2(1 + cos <j))(k - 2) - 2(q + 1) cos 0]c fc _ 2 + 
(a + 4- A:-2A)c fc _ 3 - 

The asymptotic property follows from the fact that, as in Lemma 2.1, the 
function f(x,w) can be written as f(x,w) = exp[ijj(w)], where ip(w) = 
rw 2 {aQ+a\w-\- . . .). Hence, the coefficients c k have the asymptotic behaviour 
c k = 0(rL fc / 2 J), as r — > oo. The first term approximation can be written as 
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In this case a limit for large values of r (or A or x) cannot obtained from 
the above representations. We can obtain a limit by putting A = (a + l)/2. 
Then, as follows from the recursion relation dl8|), we have Ck = 0{4> 2 ) as 
0^0, and we obtain the limit of the Askey scheme 

lim tPt +1)/2 [(a + 1)(1 - cos0) - £)/(2sin0);0] = 

q>—>0 

This includes the limit of the Askey scheme (cf . || ) 

nmP(«+ 1 )/ 2 (-e/(20);^)=LW(e). 

4.2 Two free parameters 

Next we solve c\ = 0, C2 = for A and B, with C = a. This gives (cf. flli]) ) 



A = ^/4(Acosc/> + xsin(/>) 2 - 2(a + 1)(A cos 20 + x sin 20), 

2(A cos + x sin0) + ^4 

is = , 

a + 1 

and the first term approximation can be written as 

P^(x;<P)=L^(0 + o(r n - 2 ), £ = A/B, 

as r — > oo, uniformly with respect to 0. 

As an alternative, we solve ci = 0, C2 = for ^4 and C, with B = 1. This 
gives 

A = 2[x(sin — sin 20) + A(cos — cos 20)], 

C = 2[x(2sin0 - sin 20) + A(2cos0 - cos 20)] - 1. 

and the first term approximation can be written as 

P^(x;<P)=[L^(0 + o(r n - 2 )}, £ = A, a = C. 

as r — * oo, uniformly with respect to 9. 

Solving A = £, C = a for x and A, we obtain 

A = (1 - cos0)£ + -(a + l)(2cos0 - 1), 

(18) 

2(£ - a - 1) cos 2 + (a + 1 - 2£) cos + a + 1 - £ 

x = . 

2 sin 
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Then c 3 = |(a + 1 — 2£)(1 — cos 0) and c fc = 0(0 2 ) as — ► 0, which follows 
from deriving a recursion relation for c^. 

Using these values of x and A, we obtain the limit 

limPW(x;0) = LH(O. 
4.3 Three free parameters 

We solve c\ = 0, C2 = 0, C3 = for A, B and C. This gives 

2 sin 0(x sin <f> + A cos 0) \J x 2 + A 2 2r sin sin ^ (# + 0) 



A = 



x sin 2^ + A cos 20 + sin 4>Vx 2 + A 2 sin | (0 + 30) 

x sin 20 + A cos 20 + sin 0V x 2 + A 2 _ sin | (0 + 30) 



C + l = 2 



x sin + A cos sin ^ (0 + 0) 

x sin 20 + A cos 20 + 2 sin 0Vx 2 + A 2 _ 2r [sin(0 + 20) + 2 sin 0] 



5 2 B 2 
The first coefficients Cfc are given by 

Co = 1, ci = c 2 = c 3 = 0, 

C4 = ^{sin(0 + 40) + [sin - sin(0 + 2<p)]B 2 }. 

The first term approximation can be written as 

PW(x; 0) = B» [l<? )(0 + O (r-3)] , f = | = 
as r — > 00, uniformly with respect to 9. 

5 Jacobi, Meixner and Krawtchouk to Laguerre 

We give the results for one free parameter only. 
5.1 Jacobi to Laguerre 

Let R(w) = y/1 — 2xw + w 2 . The generating function reads 

F(x,w) = i '— = \ P^' p >{x)w n . 

n n=0 
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As in Lemma 3.1, we define coefficients Ck, and the expansion reads 

n 

pM(x) = J2 B n - k c k L[%(0, ? = A/B. 

k=0 

We consider a + (3 — > oo, and solve c\ = for A, with B = 1 and C = a. 
This gives 

A=~(a + p + 2)(l-x). 
The first coefficients c k are given by 

co = l, ci = 0, c 2 = -[-a + 3/3- 2(a + 3/3 + 4)x + (3a + 3/3 + 8)x 2 ]. 

8 

The first term approximation can be written as 

A limit can be obtained by writing x = 1 — 2£/(a + (3 + 2) . Then we have 
Cj, = 0(1/ (3) as (3 — > oo for k > 2, and we obtain 

lim P^) [1 - 2£/(2 + a + /3)] = (£), 

which includes the limit of the Askey scheme (cf. ||) 
hm P^)(l-2£//3)=4*)(£). 



5.2 Meixner to Laguerre 

The generating function reads 



F(w) = ( 1 - -X (1 - «,)-^-» = £ %M n (^,c)/, 

V C / 77-. 

N ' n=0 

and we define c k as in Lemma 3.1. The expansion reads 

n 

M n (x;f3,c)=Y,B n - k c k L^ k (Z), £ = A/B. 

We solve c\ = for A, with B = 1 and C = a. This gives 

. (a-/3 + l)c+(l-c)x 
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The first coefficients c& are given by 

(1 + a - (3)c 2 + (2c - c 2 - l)x 
c = 1, ci = 0, c 2 = ^ • 

The first term approximation can be written as 

(a - /? + l)c+ (1 - c)x 



M n (x;P,c)=L( n a \t) + 0(p n - 1 ), £ 



A limit can be obtained by putting j3 = a + 1 and writing x = c£/(l — c). 
Then we have c 2 = (c — l)£/(2c), and c^ = 0(1 — c) as c — ► 1 for k > 2. We 
obtain the limit of the Askey scheme (cf. [||) 

lim M n (c£/(1 -c);a + l,c) = 



L^(0) 



5.3 Krawtchouk to Laguerre 

Let q := — -. The generating function reads 



F(w) = (1 - giu)* (1 + w) 



JV 

N-x 



2 (^W^'^ 

^0 W 



and we define c& as in Lemma 3.1. The expansion reads 



n , 



K n (x;p,N) = J2B n - k c k L^} k (0, £ = A/B. 

k=0 



We solve c\ = for ^4, with B = 1 and C = a. This gives ^4 = a + l — iV + 
(1 + g)x. The first coefficients c& are given by 

c = l, ci = 0, c 2 = ^[l + a - 3N + (3 + 2q - q 2 )x}. 
The first term approximation can be written as 

r^jK n (x;p,N) =4 c )(C)+0(iV n - 1 ) , £ = a+l-N+(l+q)x, N -> oo. 
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